Magnetic phases of one-dimensional lattices with 2 to 4 fermions per site 
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We study the spectral and magnetic properties ol one-dimensional lattices filled with 2 to 4 
fermions (with spin 1/2) per lattice site. We use a generalized Hubbard model that takes account 
all interactions on a lattice site, and solve the many-particle problem by exact diagonalization. 
We find an intriguing magnetic phase diagram which includes ferromagnetism, spin-one Heisenberg 
antiferromagnetism, and orbital antiferromagnetism. 
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I. INTRODUCTION 

Artificial lattices resemble periodic arrangements of 
quantum wells confining a small number of particles. Ex- 
perimentally, both lateral and vertical lattice structures 
can be realized. Examples are arrays of quantum dots 
in semiconductor heterostructures 1 "! 2 -^ confining the con- 
duction electrons, or optical lattices - stable periodic 
arrays of potentials created by standing waves of laser 
light 4 ^. Varying the intensity of the laser light, one 
can change the depths of the single traps, i.e. the sin- 
gle sites. In such egg-box like potentials, experimental- 
ists can confine ultra-cold atoms, of bosonic or fermionic 
character^&^ i 10 : 11 , achieving particle numbers on the 
sites that are even less than three. The strengths and 
even the sign of the interactions between the atoms can 
be tuned by Feshbach resonance a 12 ' 13 i 14 ' 15 i 16 ' 17 . 

The basic difference between artificial lattices and nor- 
mal lattices (such as the crystal structure of solids) is, 
that in artificial lattices the particles confined in the lat- 
tice do not play any role for determining the intrinsic 
lattice structure. A possible degeneracy of the many- 
particle states can then not be removed by lattice distor- 
tion. Instead, it may lead to internal symmetry breaking 
and, for example, to spontaneous magnetism and super- 
conductivity. Recent experiments have inspired much 
theoretical work on artificial lattices, both with cold 
atomsi&i&^i 2 ^ and quantum dots 2 ^!. 

Mean-field calculations based on the spin-density func- 
tional theory predict that Hund's first rule determines 
the total spin of an isolated, individual lattice sit o 25 ' 26 . 
The magnetism of the lattice as a whole then depends on 
the total spin of the individual lattice sites, on the lattice 
structure and on the coupling between the site s 27 ' 28 : 29 ' 30 . 
A simple tight-binding model with a few parameters can 
account for most of the these findings 3 -!. Related results 
have been obtained for quantum dot molecules using the 
density functional method 32 . 

The eigenstates of single quantum dots with a few elec- 
trons can be calculated "exactly" (i.e. to a high degree of 
convergence with respect to the necessary restrictions in 
Hilbert space) by diagonalizing the many-body Hamilto- 
nian (for a review see Ref.2 6 -). Methods beyond the mean- 



field approximation have also been applied to quantum 
dot molecules^i^ 3 ^ 7 -^. 

For a lattice with strongly correlated particles the 
generic model is the Hubbard model, which has been 
amply studied in the case of one state per lattice site 
(for reviews se o 39 ' 40 ). From an experimental viewpoint, 
it has been argued that the Hubbard approach is ideal 
for describing contact-interacting atoms in an optical lat- 
tice^ ' 41 ' 42 ' 43 ' 44 . The one-dimensional Hubbard model is 
exactly solvable using the Bethe ansata 4 ^.. The mag- 
netism of finite molecule s 46 ' 47 and quantum rings 4 - have 
also been studied in the simple Hubbard model. 

The purpose of this paper is to study the magnetism 
of an artificial one-dimensional (ID) lattice in the case 
where the lattice site is filled on average with 2 to 4 
fermions, which can be electrons in a quantum dot lat- 
tice or fermionic atoms in an optical lattice. We call the 
electrons or atoms generally as particles. We assume the 
confining potential in each lattice site to be quasi-two- 
dimensional and nearly harmonic at the bottom. In this 
case the ls-level of each lattice site is filled, and the de- 
generate lp level is partially filled. We use a generaliza- 
tion of the Hubbard model to describe the interactions: 
The particles interact only within a lattice site. We solve 
the Hubbard Hamiltonian by exact diagonalization for a 
finite length of the lattice using periodic boundary con- 
ditions. The results show many different magnetic struc- 
tures which are analyzed through their relations to the 
Heisenberg model and the simple single-state Hubbard 
model. 



II. THEORETICAL MODELS 
A. ID lattice with p-orbitals 

We consider an artificial lattice where the confining po- 
tential at each lattice site is nearly harmonic and quasi- 
two-dimensional so that the single-particle level structure 
in each lattice site is Is, lp, 2s\d etc. We assume that in 
all cases the ls-state is filled completely, and the doubly 
degenerate lp-level is partially filled. Furthermore, we 
assume that the shells are well separated in energy, so 
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FIG. 1: Schematic pictures of ID lattices considered. Each 
lattice site has p x and p y orbitals. In the lateral case these 
are shown as light and dark-gray densities. Here, the hop- 
ping probabilities t and ti between neighboring lattice sites 
are different for p x and p y orbitals. In the vertical case it is 
natural to use states with 'rotating orbitals' p+i and p_i with 
circularly symmetric densities. In this case there is only one 
hopping probability t. 



that we can neglect the mixing of the lp shell with the 
Is or 2sld shells. This leads to a generalized Hubbard 
model which has in each lattice site only two orbitals 
which we call either p x and p y or p_i and p+i, respec- 
tively. The latter notation refers to orbitals with angular 
momentum quantum numbers —1 and +1 (clockwise or 
counterclockwise rotation of the p-state) . 

The two kinds of ID lattices considered are schemat- 
ically shown in Fig. [TJ In the case of semiconductor 
quantum dots, these are often called lateral and vertical 
structures. In the lateral lattice the hoppings between 
neighboring p x and p y states are different and denoted 
by t and t2, where t% < t. For the vertical lattice, it 
is natural to use the angular momentum states p-\ and 
p+i- In that case there is only one hopping parame- 
ter t (or equivalently t^ = i). Note that for the single- 
particle wave functions we have "0+1 = (V-Vx + * - 0p B )/v / 2 
and ip-i = (ip Px - itp Py )/y/2. 



B. Hubbard model 

We assume a generalized Hubbard model Hamiltonian 

H = J + U, (1) 

where the first term represents inter-site hoppings be- 
tween neighboring lattice sites and the second term intra- 
site two-body interactions. 

Hoppings preserve spin, and are equal for spin-up and 
spin-down particles. Thus, J separates into two symmet- 
ric spin parts: J = 53cr=T 1 CT ' ^ or our one ~dimensional 
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where n is the lattice site index, and j and j' denote 
the p-orbital in question. (In the simple Hubbard model, 
there would be only one space state per site, and the 
j-indices not needed). 

Some of the hopping integrals Jjy are zero due to sym- 
metry. The non-zero integrals are treated as essentially 
free model parameters, t and ti- Thus, we have 
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for the lateral and vertical lattice, respectively. (Note, 
that in the case £2 = t, the lateral model actually is iden- 
tical to the vertical model, irrespective of the different 
p-orbit basis used). 

We approximate the two-body interactions in the spirit 
of the tight-binding model: The particles only interact 
when they are at the same lattice site. Thus, U separates 
in the symmetric parts representing interactions on each 
site n: U = J2 n &n- Within a site, full (spin-indpendent) 
two-body interaction is allowed, which yields 
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where Uj 1 j 2 j 3 j i are the direct space matrix elements of 
on-site interaction, depending on the interaction itself 
and the j-orbits in question, i.e. the eigenstates of the 
confining potential. 

For contact interactions, the ratios of the different ma- 
trix elements are independent of the confining potential, 
as long as it has circular symmetry. For the non-zero ma- 
trix elements (together with those obtained by allowed 
j-index permutations), we obtain 



3i 


h 


33 


H 




Px 


Px 


Px 


Px 


3/7 


Px 


Px 


Py 


Py 


U 


Px 


Py 


Px 


Py 


U 


Px 


Py 


Py 


Px 


u 


Py 


Py 


Py 


Py 


3U 


P-l 


P-i 


P-i 


P-i 


2U 


P-l 


P+i 


P-i 


P+i 


2U 


P-l 


P+i 


P+i 


P-i 


2U — A 


P+l 


P+i 


P+i 


P+i 


2U 



where U is the only parameter describing the strength 
of the interaction. All together, we thus have three pa- 
rameters t, ti and U . One of them can be fixed to set 
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the energy scale. We choose this to be t and represent 
the results for t = 1 (all energies are given in units of t). 
In some cases with vertical lattices we also consider an 
interaction of finite width. This can be mimicked by de- 
creasing one of the matrix elements by a small amount A, 
as indicated in the above table. For contact interactions, 
A = 0. 

We solve the Hamiltonian for a lattice with L lattice 
sites using periodic boundary conditions (J connects also 
the last and the first site). The Lanczos method is used 
to find the low-energy eigenvalues and eigenvectors of the 
Hamiltonian matrix. We take advantage of the periodic- 
ity of the lattice and solve the eigenvalues separately for 
each Bloch fc-value. In practice this means that, instead 
of using "site-states" | nja) as a single-particle basis to 
span the Fock space, one uses Bloch states of the tight- 
binding model (eigenstates of J). In this study, the hop- 
ping does not mix the p x and p y orbitals in the lateral 
case, nor p_\ and p +1 orbitals in the vertical case. We 
then have separate, simple bands with energy eigenvalues 

et (*) = -2tco S (^V (4) 

where k takes integer values 0, 1, • • • , L — 1. Note that 
in the lateral case, the p x and p y bands have different 
widths for t and t%, respectively. In the vertical case, the 
widths are always the same. 

We do not take advantage of the fact that the Hamil- 
tonian does not depend on spin, but diagonalize the sys- 
tem for S z = and only afterwards determine the total 
spin S for each many-particle state. The total number 
of particles is denoted by N and the numbers of spin- 
up particles and spin-down particles by iVj and N^. We 
note that because of the spin degree of freedom, the max- 
imum number of particles in a lattice with length L, is 
^Vmax = 4L. The filling fraction by v = N/L gets values 
from to 4. In this study we consider only the region 
v = • • • 2. Due to the symmetry of the Hamiltonian the 
region v = 2 • • • 4 will have similar properties. 

As discussed earlier, we assume that the hopping can 
occur only between the nearest neighbours. It should be 
noted, however, that the interaction part of the Hamilto- 
nian allows intra-site hopping, via scattering from one 
single-particle state to another inside any lattice site. 
This becomes important especially in the case where 
t-2 = 0, where the hopping only occurs through the p x 
states. 



C. Heisenberg model 

It is well-known that the simple Hubbard model in 
the limit of large U/t approaches the antifcrromagnetic 
Heisenberg model. In this case, the low-lying eigenstates 
are characterized by one spin 1/2 particle on each site. 
In a similar way, in some limiting cases, our results with 
p-orbitals approach those of the Heisenberg model with 
S = 1 (two particles on each site with aligned spins), 



or with 5=1/2 (polarized system with one fermion on 
each site, with the p-orbitals playing the role of the spin 
components). The effective Hamiltonian is then 

H c s = Jeff ^ S„ ■ S n+ i + constant, (5) 

where J c g is the effective exchange interaction and S« 
the spin operator for site n. We compare the Heisenberg 
and Hubbard model for the case of four sites, L = 4, 
where the spectrum of the antiferromagnetic Heisenberg 
model is exactly solvabl o 48 ' 49 . 

III. RESULTS 

A. A single lattice site with two particles 

A single site with two particles obeys Hund's first rule 
to maximize the spin. The energy difference between the 
lowest S = and S = 1 states is the 'exchange splitting' 
and equals AE = Es=o — Es—i = 2U. In the case of 
a finite-range interaction the exchange splitting is 2U — 
A. Table U gives the energy spectrum of a single lattice 
site. We will see below that in the limit of large U the 
half-filled system (y = N/L — 2) becomes a Heisenberg 
antifcrromagnet with 5 = 1. 

TABLE I: Energy levels and corresponding total spin of a 
single lattice site with two particles. 



State No 


E 


3 


1 


A 


1 


2 


2U 





3 


2U 





4 


4(7 - A 






B. Half-filled vertical lattice: N — 2L 

In the half-filled case there is one particle per orbital. 
When U is large, each lattice site will have spin S = 1 due 
to the large exchange splitting. The only way to allow 
particles to hop from one site to the neighboring one is to 
orient the total spins of neighboring sites opposite, i.e. 
with antiferromagnetic order. For ferromagnetic order, 
the hopping would be prohibited by the Pauli exclusion 
principle. In this case the total energy of the system 
would be zero (assuming A = 0). For antiferromagnetic 
order the allowed hopping can reduce the energy to a 
slightly negative value. 

We will first study the case of vertical lattice (t^ = t) . 
Figure [2] shows the low-energy energy levels for four lat- 
tice sites and eight particles (L = 4, JV = 8), calculated 
for different values of U. All the levels with energy E~ U 
are shown. For the largest value U = 50 the spectrum 
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FIG. 2: Low-energy levels (_B < U) for L = 4 and N = 8 
calculated for different values of (7. The numbers next to 
the levels denote the total spin S of the many-particle state. 
The wave vector k has values 0, 1, and 2. The symbols plus, 
square, star, circle and dot correspond to U values 2, 2.5, 
5, 10, and 50, respectively. The energy levels for U = 50 
agree with those of the S — 1 Heisenberg model with 0.01 % 
accuracy. 

agrees with that of the Heisenberg model (J e g = 2/J7) 
with 0.01 % accuracy. The Heisenberg model for four 
sites is an exacly solvable textbook problem 48,49 . It is 
interesting to notice that even for U = 2 the spectrum 
is qualitatively still the same. Only when U < 1.5 new 
states start to appear in the low-energy spectrum. 

C. Vertical lattice polarized fermions: The 
noninteracting case 

We will now consider polarized fermions (e.g. elec- 
trons or fermionic atoms with TV = N-\). For contact 
interactions between the fermions, the problem becomes 
non-interacting since the Pauli exclusion principle forbids 
two fermions to be at the same state. The energy spec- 
trum can then be constructed by filling particles to the 
Bloch states (Eq. (QJ) which are solutions of the nonin- 
teracting Hamiltonian J. 

In this (trivial) case, it is important to note that (i) 
each single-particle Bloch state is doubly degenerate due 
to the two states per site, and (u) only for particle num- 
bers N = An + 2 the ground state is non-degenerate (n 
is a non- negative integer), (ii) implies that the ground 
state energy (of polarized fermions) as a function of N 
has local minima for N = 2, 6, 10, • • • . We will see later 
that these special values form single domain ferromagnets 
when N <2L. 



D. Orbital antiferromagnet of polarized fermions: 
Vertical lattice with N = L 

Let us now consider polarized fermions with a finite 
range interaction and one fermion per lattice site (N — 
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FIG. 3: Low-energy levels (E < 0) for polarized fermions 
with L = 4 and iV = 4 calculated for different values of A 
for U = 50. The numbers next to the levels denote the total 
spin of the corresponding state of the spin- 1/2 Heisenberg 
model. The wave vector k has values 0, 1, and 2. The symbols 
star, circle and dot correspond to A values 2, 10, and 50, 
respectively. The energy levels for A = 50 agrees with those 
of the 5=1/2 Heisenberg model with 0.1 % accuracy. 



AT-f = L). The finite range here means only that A > 0. 
However, the finite range does not lead to interaction 
of particles sitting at different lattice sites. Each site 
still has two p states. The large U limit in this case 
is an antiferromagnet where the 'magnetic moment' in 
each lattice site is not the spin but the orbital angular 
momentum of the p states, which can have the two values 
+1 or -1. 

There are two reasons for this state to become the 
ground state. First, it costs energy (by A) for two par- 
ticles to occupy the same site. Thus, the particles prefer 
to be at different sites. Second, the particles can only 
hop to the neighboring site if they are at different orbital 
states. Although the particles prefer to be at different 
sites, a small amount of 'virtual' hopping is necessary to 
reduce the energy. 

Again, we compare the spectrum with the exact result 
of the Heisenberg model for four particles. In Fig. [3l all 
the low energy states are plotted for different k values. 
For large U and A the agreement between the Hubbard 
model and Heisenberg model becomes perfect with J e g = 
1/A. 



E. Vertical lattices with iV < 2L: Ferromagnetism 

Next, we consider vertical lattices with contact inter- 
actions and large values of U (U > 10). The results show 
that the ground states for N — 2, 6, and 10 have max- 
imum possible total spin, i.e. they are ferromagnetic. 
This is true for all values of L > N/2 where the com- 
putations could be performed (the matrix sizes increase 
very fast with L). The ferromagnetic ground state can 
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be understood as follows. When L 3> N/2 the ferromag- 
netic state allows particles to move freely in the lattice, 
as even in cases where two particles are in the same lat- 
tice site, they do not interact. In other words, particles 
with the same spin can pass each other without any cost 
of energy. If the particles have opposite spin, however, 
they suffer repulsive interaction whenever they are at the 
same lattice site - even if they are at different p-states. 
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FIG. 4: Motion of holes in the nearly half-filled case with large 
U which prevents two opposite spins to be at the same lattice 
site. In the ferromagnetic case, the holes can move indepen- 
dently, while in the antiferromagnetic case they are bound 
together since their separation costs energy, as indicated by 
U. 

Let us now consider what happens if we start from 
the antiferromagnetic L = N/2 state and increase L by 
one. Also in this case the results show that, indepen- 
dent of the particle number (here, 2 < N < 10), the 
ground state is ferromagnetic. In the antiferromagnetic 
case with L = N/2 the lowest energy is proportional to 
— 1/U, which for large U is very small. The energy of the 
ferromagnetic case with L = N/2 is zero since there is no 
room for hopping (the single-particle bands are filled) . If 
now one lattice site is added, the ferromagnetic energy 
becomes —At. This is because there are now two freely 
moving holes in the system, as illustrated in Fig. SJ The 
situation is different if the system remains antiferromag- 
netic. Also in this case there are two holes, but now they 
are bound together, since their separation costs energy, 
as illustrated in Fig. 2J The total energy of the antiferro- 
magnetic state will necessary be above the ferromagnetic 
energy —At. Consequently, adding one lattice site to the 
antiferromagnetic L = N/2 case transforms it to a fer- 
romagnetic state. Alternatively, we can start from the 
half-filled case and remove one particle. In the ferromag- 
netic case the hole is free and has an energy of —2t, while 
in the antiferromagnetic case the hole is localized and its 
energy is zero. 

As mentioned above, the ferromagnetic ground state 
has total spin S = N/2 for N = 2, 6, 10, • • • . However, 
the situation is more complicated for N = A, 8, 12, 
In these cases the total spin of the ground state is S = 
for all L > N/2. Nevertheless, we argue that also these 
cases are ferromagnetic, but now the ground state has a 
spin wave which rotates the spin once within the length 
L. Alternatively, we can apply the picture that the fer- 
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FIG. 5: Pair-correlation of vertical lattices of different lengths 
L — 4, 6, 8 with N — 8 fermions. In each case the reference 
site is where one particle with spin-up is fixed. The solid 
line with black dots shows the total up-up correlation and the 
dashed line with open circles the up-down correlation. The 
uppermost panel shows the antiferromagnetic ground state of 
L — 4, the center panel and the lowest panel show the two 
domains of the ferromagnetic ground states of L = 6 and 
L = 8, respectively. 



romagnetic ground state consists of two domains with 
opposite spin directions. The reason for this behavior 
is the fact that for these particle numbers the ferromag- 
netic state is degenerate and the spin wave (or domain 
formation) provides a way to remove the degeneracy and 
reduce the total energy. 

Figure [5] shows the pair-correlation function of N = 8 
particles for L = 4, 6, and 8. We fix one particle in a 
state, say p+i with spin-up in lattice site and determine 
the conditional propability of finding the other spin-up 
and spin-down particles on the other lattice sites. Figure 
[5]shows clearly that for L = 4 = N/2 the result is antifer- 
romagnetic, while for L — 6 and L — 8 the spin changes 
direction only once within the length L, as it would hap- 
pen for the longest possible spin wave. It is interesting 
to note that in fact, the system with two states per site 
is very different from that with only one s state per site. 
In the latter case the system remains antiferromagnetic 
(for large U) for all values of Zi 48 i 50 . 
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F. Lateral lattices: t^ < t 

In the lateral lattice, as shown in Fig. [T] the hopping 
parameters t and t 2 for the two p-states are different. 
The structure of the ground state and the many-particle 
spectrum then depends on the ratio t 2 /t. We will now 
study the magnetism as a function of this ratio and of 
the filling fraction N/L. 

For different values of t and t 2 , for noninteracting par- 
ticles we have two cosine bands, Eq. ((4]), which reach 
from — 2t to +2t and from — 2t 2 to +2t2, respectively. 
Lets consider first the ferromagnetic case with low filling 
and remember that for contact interactions the system 
becomes non-interacting. In the limit of low filling and 
t 2 < t, only the t-band is occupied. This is equivalent to 
the simple one-state Hubbard model. But we know that 
the ground state of the one-state Hubbard model is an- 
tiferromagnetic in the case of low filling. Consequently, 
the ground state will be antiferromagnetic whenever the 
corresponding ferromagnetic state would only occupy the 
i-band. This condition can be easily derived, 

N 1 , , . . 

A similar argument can be used to show that for 

N/L> 2- J- cos" 1 (ta/t) (7) 

Z7T 

the system is also antiferromagnetic. Here, the holes in 
the ferromagnetic case only occupy the i-band. In be- 
tween these limits, both bands are partially filled and 
determining the magnetism is more complicated. For 
t 2 = t the lateral lattice equals the vertical lattice. We 
have shown above that this case should always be ferro- 
magnetic. We can thus expect that for t 2 close to t, the 
system is ferromagnetic. 

Figure[5]shows the magnetism of the ground state, as a 
function of both the number of particles per site (on the 
p-states) and the ratio t-ijt of the two hopping param- 
eters, calculated by diagonalizing the Hubbard Hamil- 
tonian. The figure also shows the limits given by Eqs. 
^ and ((7]). Indeed, we see that between these lim- 
its, the ground state is mainly ferromagnetic, while out- 
side these limits it is always antiferromagnetic. Figure [5] 
shows results computed for 2, 6 and 10 particles, where 
the ferromagnetic phase is simple and seen as the spin 
being at maximum S = N/2. As discussed above, for 
N = 4, 8, 12, ■ • • the ferromagnetic state has a spin- 
wave (or domains) and interpretation of the magnetic 
structure is more difficult. Nevertheless, results com- 
puted for those particle numbers seem to agree with the 
phase diagram shown in Fig. [5J The results in Fig. [5] are 
computed for U = 10. We repeated some of the points for 
larger values of U and found the same magnetic states. 

It is interesting to compare the above results of the gen- 
eralised Hubbard model with those of the density func- 
tional mean field theor y 28 i 29 i 30 . The qualitative agree- 
ment is perfect: In the case of two p-particles per site 
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FIG. 6: Magnetic phase diagram of the one-dimensional lat- 
eral lattice. The vertical axis shows the number of particles 
per site (on p-states) and the lateral axis the ratio ti/t of the 
two hopping parameters. The dashed lines show the borders 
between which the narrower band {t% band) is filled. Outside 
this area the lattice is antiferromagnetic (except at t-z = t). 
The filled (open) symbols show the ferromagnetic (antiferro- 
magnetic) ground states obtained with exact diagonalization 
of the Hubbard Hamiltonian for 10 (triangles), 6 (squares), 
and 2 (circles) particles. Crosses show those results for N = 6 
which are not ferro or antiferromagnetic, i.e. < S < N/2. 
The numerical results are for U = 10. 



(N/L = 2) the system shows antiferromagnetism of spin- 
one quasiparticles, while in the case of one p particle per 
site (N/L = 1) the system is ferromagnetic. An even 
simpler tight-binding model 27 also gives a similar phase 
diagram. Due to the symmetry of the Hamiltonian, it is 
natural that also above the filling N/L — 2 one obtains 
a ferromagnetic region with its center at N/L = 3. 

For small values of t 2 the corresponding single-particle 
band becomes very narrow. In this case the ferromag- 
netism can be understood with the Stoner mechnanism 51 : 
The Fermi level is in the region of large density of states 
and induces a ferromagnetic state. In ID this effect is 
particularly strong due to the singularities in the density 
of states 2 ^. 



IV. CONCLUSIONS 

We studied the magnetism of one-dimensional artificial 
lattices made of quasi- two-dimensional potential wells, 
for up to four particles per lattice site, i.e. in the region 
where the lp level is filled. We froze the Is particles 
and considered only the lp states. Numerical diagonal- 
ization of a generalized Hubbard model was performed 
for several particle numbers and filling fractions. The re- 
sults were analyzed using the antiferromagnetic Heisen- 
berg model and single-particle models. 

In the resulting phase diagram, the vertical lattice is 
ferromagnetic, except at a singular point with exactly two 
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p-type particles per site. For lateral lattices the ground 
state is antiferromagnetic for small fillings and close to 
half-filling of the p-shell, but ferromagnetic around the 
region with one p-particle per site. A simple model for 
the ferromagnetic region was suggested. 

If the particle number is a multiple of four (N = 
4,8, 12, •••), the ferromagnetic state has a spin wave 
which removes the degeneracy and yields a total spin 
S = Q. 



For polarized fermions the half-filled case shows "or- 
bital antiferromagnetism" where in successive lattice sites 
the particles rotate clockwise and counter-clockwise. 
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